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DESCRIPTION 



STOCHASTIC CODEBOOK EXCITATION VECTOR CODING METHOD 



5 Technical Field 

The present invention relates to a stochastic 
codebook excitation vector coding method in a CELP speech 
coding apparatus/speech decoding apparatus. 



10 Background Art 

When speech signals are transmitted in a packet 
communication system typified by Internet communication^ 
a mobile communication system, or the like, compression 
and coding techniques are used to improve the speech signal 

15 transmission efficiency. Many speech coding methods 
have been developed to date, and many low bit rate speech 
coding methods developed in recent years, such as CELP, 
separate a speech signal into spectrum envelope 
information and spectrum detailed structure information, 

20 and perform compression and coding of the separated 
information. 

In a CELP speech coding apparatus, synthetic speech 
vectors are calculated for all combinations of adaptive 
code vectors stored by an adaptive codebook and fixed 
25 code vectors stored by a stochastic codebook, distance 
calculation is performed for each synthetic speech and 
input speech signal, and the adaptive code vector index 
and fixed code vector index for which the distance is 
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smallest are found. 

One known stochastic codebook is an algebraic 
codebook. This codebook enables a stochastic codebook 
search to be performed with a comparatively small amount 
5 of calculation, and has consequently been widely used 
in CELP in recent years. 

An excitation vector of an algebraic codebook is 
composed of a small number of pulses with an amplitude 
of 1 and polarities ( + , -) , and the pulses (in this case, 
10 excitation vector waveform candidates) are positioned 
so as not to overlap each other. 

For example, when the subframe length is 32 and the 
number of pulses (= number of channels) is 4, the number 
of pulses per channel is 32/4 = 8, and the channel 0 pulse 
15 positions iciO[iO], channel 1 pulse positions icil[il], 
channel 2 pulse positions ici2[i2], and channel 3 pulse 
positions ici3Ii3] are as shown below. Here, iO, il, 12, 
and 13 denote indexes of the respective channels. 

iciO[iO] = {0, 4, 8, 12, 16, 20, 24, 28} 
20 icil[il] = {1, 5, 9, 13, 17, 21, 25, 29} 

ici2[i2] = {2, 6, 10, 14, 18, 22, 26, 30} 

ici3[i3] = {3, 7, 11, 15, 19, 23, 27, 31} 

A conventional stochastic codebook codes the pulse 
positions of each channel independently, and takes codes 
25 combining these with polarity codes as stochastic 
excitation vector codes. 

For example, in the above case of a subframe length 
of 32 and 4 channels, a convent ional codebook 103 
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represents a pulse position of each channel as 3 bits, 
and together with the polarity code, performs coding using 
a code of (3+l)x4=16 bits. 

However, a problem with the above conventional 
5 stochastic codebook coding method is that, if the bit 
rate is low the bits assigned to each channel are also 
limited, and there are positions where there is no pulse 
at all, so that variations of an excitation vector waveform 
corresponding to a code (position information) decrease, 
10 and sound quality degradation occurs. 

In the above case of a subframe length of 32 and 
4 channels, for example, there are positions where there 
is no pulse at all if coding is performed with fewer than 
16 bits. 

15 

Disclosure of Invention 

It is an object of the present invention to provide 

a stochastic codebook excitation vector coding method 

that enables variations to be secured so that there are 
20 no positions where there is no pulse at all while achieving 

a reduction of the number of bits used when coding 

stochastic codebook pulses. 

This object is achieved by associating a pulse 

position of a predetermined channel with a pulse position 
25 of another channel, searching for a pulse position by 

means of a predetermined algorithm, and taking a found 

pulse position code and a polarity code as a stochastic 

excitation vector code. 
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Brief Description of Drawings 

FIG.l is a block diagram showing the configuration 
of a CELP speech coding apparatus; 
5 FIG. 2 is a flowchart showing an example of a pulse 

search algorithm for each channel in a coding method 
according to Embodiment 1 of the present invention; 

FIG. 3 is a flowchart showing an example of a pulse 
search algorithm for each channel in a coding method 
10 according to Embodiment 1 of the present invention; 

FIG. 4 is a flowchart showing an example of a pulse 
search algorithm for each channel in a coding method 
according to Embodiment 2 of the present invention; and 
FIG. 5 is a flowchart showing an example of a pulse 
15 search algorithm for each channel in a coding method 
according to Embodiment 2 of the present invention. 

Best Mode for Carrying out the Invention 

FIG.l is a block diagram showing the configuration 

20 of a CELP speech coding apparatus . An input speech signal 
is input sequentially to the speech coding apparatus 
divided into processing frames at time intervals of 
approximately 20 ms . 

The input speech signal input to the speech coding 

25 apparatus every processing frame is first supplied to 
an LPC analysis section 101. LPC analysis section 101 
performs LPC (Linear Predictive Coding) of the input 
speech signal and obtains an LPC coefficient, performs 

L 
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vector quantization of the LPC coefficient to produce 
an LPC code, and decodes this LPC code to obtain a decoded 
LPC coefficient. 

An excitation vector creation section 104 reads an 
5 adaptive code vector and fixed code vector respectively 
from an adaptive codebook 102 and stochastic codebook 
103, and sends these to an LPC combining section 105, 
LPC combining section 105 performs combining filtering 
of the adaptive code vector and fixed code vector supplied 

10 from excitation vector creation section 104, and the 
decoded LPC coefficient provided from LPC analysis 
section 101, with an all pole type combining filter in 
the filter coefficient, and obtains a combined adaptive 
code vector and combined fixed code vector. 

15 A comparison section 106 analyzes the relationship 

between the combined adaptive code vector and combined 
fixed code vector output from LPC combining section 105, 
and finds adaptive codebook optimum gain to be multiplied 
by the combined adaptive code vector, and stochastic 

20 codebook optimum gain to be multiplied by the combined 
fixed code vector. 

Comparison section 106 also adds together the vector 
obtained by multiplying the combined adaptive code vector 
by the adaptive codebook optimum gain and the vector 

25 obtained by multiplying the combined fixed code vector 
by the stochastic codebook optimum gain, and obtains a 
combined speech vector, and performs a distance 
calculation on the combined speech and input speech signal . 
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Then comparison section 106 obtains the adaptive code 
vector stored by adaptive codebook 102 and the combined 
speech vector stored by stochastic codebook 103, and finds 
the adaptive code vector index and fixed code vector index 
5 for which the distance between the combined speech and 
input speech signal is smallest. Comparison section 106 
then sends the indexes of the code vectors output from 
the codebooks, the code vectors corresponding to the 
respective indexes, and the adaptive codebook optimum 

10 gain and stochastic codebook optimum gain, to a parameter 
coding section 107. 

Parameter coding section 107 codes the adaptive 
codebook optimum gain and stochastic codebook optimum 
gain and obtains a gain code, and outputs the gain code, 

15 the LPC coefficient provided by LPC analysis section 101, 
and the indexes of each codebook together for each 
processing frame. 

Parameter coding section 107 also adds together the 
two vectors comprising the vector obtained by multiplying 

20 the adaptive code vector corresponding to the adaptive 
codebook index by the adaptive codebook gain 
corresponding to the gain code, and the vector obtained 
by multiplying the fixed code vector corresponding to 
the stochastic codebook index by the stochastic codebook 

25 gain corresponding to the gain code, and obtains a drive 
excitation vector, and updates the old adaptive code 
vector in adaptive codebook 102 with the drive excitation 
vector . 
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Combining filtering by LPC combining section 105 
generally makes combined use of a linear predictive 
coefficient^ a high emphasis filter, and a weighting 
filter that uses a long-term predictive coefficient 
5 obtained by long-term predictive analysis of input 
speech . 

Adaptive codebook and stochastic codebook optimum 
index searches, optimum gain calculation, and optimum 
gain coding processing are generally carried out in 

10 sub frame units resulting from further division of a frame . 

In a speech decoding apparatus (decoder) , the same 
configuration of LPC analysis section 101, adaptive 
codebook 102, stochastic codebook 103, excitation vector 
creation section 104, and LPC combining section 105 is 

15 provided as shown in FIG.l, and an excitation vector 
waveform is obtained by decoding codes transmitted from 
a speech coding apparatus. 

In order to reduce the amount of calculation, 
comparison section 106 usually searches for an adaptive 

20 codebook 102 excitation vector and stochastic codebook 
103 excitation vector by means of an open- loop procedure. 
This open-loop search procedure is described below. 

(1) First, excitation vector creation section 104 
chooses excitation vector candidates (adaptive 

25 excitation vectors) in succession from adaptive codebook 
102 only, LPC combining section 105 creates a composite 
tone, and comparison section 106 carries out a comparison 
of the input speech and composite tone and selects the 
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optimum adaptive codebook 102 code. At this time, gain 
is selected on the assumption that it is the value at 
which coding distortion is minimal (optimum gain) . 

(2) Next, the above-described adaptive codebook 
5 code is fixed, excitation vector creation section 104 
successively selects the same excitation vector from 
adaptive codebook 102 and stochastic codebook 103 
successively selects the excitation vector (stochastic 
excitation vector) corresponding to the comparison 

10 section 106 code, LPC combining section 105 generates 
composite tones, and comparison section 106 compares the 
sum of both composite tones with the input speech and 
determines the optimum stochastic codebook 103 code. As 
in (1) above, gain is selected at this time on the 

15 assumption that it is the value at which coding distortion 
is minimal (optimum gain) . 

Use of the above procedure to search for the optimum 
excitation vector results in a slight degradation of 
coding capability, but also a maj or reduction in the amount 

20 of calculation, compared with the method of searching 
for the optimum excitation vector by comparing 
combinations of all excitation vectors or both codebooks . 

The stochastic codebook 103 excitation vector 
search method will now be described in detail. 

25 Excitation vector code derivation is carried out 

by searching for the excitation vector that minimizes 
coding distortion E in Equation (1) below. In Equation 
(1), X denotes the coding target; p, adaptive excitation 
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vector gain; a weighting combining filter; a, an 

adaptive excitation vector; q, stochastic excitation 
vector gain; and s, a stochastic excitation vector. 
E = |x - (pHa+qHs)|2 
5 ... Equation ( 1 ) 

As the adaptive excitation vector search is 
performed by means of an open-loop procedure, stochastic 
codebook 103 code derivation is performed by searching 
for the excitation vector that minimizes coding 
10 distortion E in Equations (2) below. In Equations (2), 
y denotes the stochastic excitation vector search target 
vector . 

y = X - p H a 

E = I y - q H s I 2 

. . . Equations ( 2 ) 
15 Here, gain values p and q are determined after the 

excitation vector search, and by making gain p = gain 
q = 1, Equations (2) above can be written as Equations 
(3) below. 

X • H a 

y = X H a 

I H a M 

y H s 

E = I y H s I 2 

I H s M 

20 ... Equations (3) 

Minimizing this distortion expression is equivalent 
to maximizing function C in Equation (4) below. 
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( y H • s ) 2 

C = 

8 H H s 

. . . Equation ( 4 ) 
Therefore, in the case of a search for an excitation 
vector composed of a small number of pulses such as an 
5 algebraic codebook excitation vector, calculating yH and 
HH beforehand enables function C above to be found with 
a small amount of calculation. 

yH can be found by reversing the order of vector 
y and convoluting matrix H, and then reversing the order 
10 of the result, and HH can be found by multiplication of 
the matrices . 

Stochastic codebook 103 searches for and codes a 
stochastic excitation vector using the procedure 
described in (1) through (4) below. 
15 (1) First, as preliminary processing, vector yH and 

matrix HH are found. 

(2) Next, pulse polarities are determined from the 
polarities (+ -) of vector yH elements. Specifically, 
the polarity of the pulse at each position is matched 

20 to the value of that position in yH, and the polarity 
of the yH value is stored in another array. After the 
polarities of all positions have been stored in another 
array, yH values are all made absolute values and converted 
to positive values. HH values are also converted in 

25 accordance with these polarities by performing polarity 
multiplication. 

(3) Next, function C shown in Equation (4) is found 
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by adding yH and HH values using an n-fold loop (where 
n is the number of channels), and the pulse positions 
of the channels at which this value is largest are found. 

( 4 ) The found pulse position of each channel is coded, 
and a code combining this with a polarity code is taken 
as the stochastic excitation vector code. 

With reference now to the accompanying drawings, 
stochastic codebook excitation vector coding methods 
according to embodiments of the present invention will 
be explained in detail below. In the descriptions of 
these embodiments , an algebraic codebook is used for which 
the subframe length is 32 and the number of pulses (= 
number of channels) is 4. 

15 (Embodiment 1) 

In Embodiment 1, a case is described in which an 
index of a predetermined channel is changed in accordance 
with another channel. 

In this embodiment, channel 0 pulse positions 
20 iciO[iO], channel 1 pulse positions icil[jl], channel 
2 pulse positions ici2 [ j2] , and channel 3 pulse positions 
ici3[j3] are as shown below. 

iciO[iO] = {0, 4, 8, 12, 16, 20, 24, 28} 
icil[jl] = {1, 5, 9, 13, 17, 21, 25, 29} 
25 ici2[j2] = {2, 6, 10, 14, 18, 22, 26, 30} 

ici3[j3] = {3, 7, 11, 15, 19, 23, 27, 31} 
Here, iO (0<i0<7) is the index of channel 0, jl 
(0<jl<7) is the index of channel 1, j2 (0<j2<7) is the 
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index of channel 2, and j3 (0^j3^7) is the index of channel 
3. 

For example, the iO=0 pulse position is {0}, the 
iO = l pulse position is {4}, and so on; and the jl = 0 pulse 
5 position is {1}, the jl = l pulse position is {5}, and so 
on . 

Channel 1, channel 2, and channel 3 pulses are 
grouped into pairs. For example, for channel 1, pulses 
are grouped into group 0 {1, 5}, group 1 {9, 13}, group 
10 2 {17, 21}, and group 3 {25, 29}. 

Then, if il (0^11^3) is designated the channel 1 
group index, 12 (0^i2<3) is des ignated the channel 2 group 
index, and 13 (0^13^3) is designated the channel 3 group 
index, the relationship between indexes jl, j2, and j3 
15 and group indexes il, 12, and 13 is as shown in Equations 
( 5 ) below . 

j 1 = i 1 X 2 + ( i 0 % 2 ) 
j2=i2x2+((iO+i1)%2) 
j 3 = i 3 X 2 + (( i 1 + i 2 ) % 2 ) 
20 ... Equation ( 5 ) 

In Equations { 5 ) , the ^'%" symbol denotes an operation 
that finds the remainder when the numeric value on the 
left of (index) is divided by the numeric value on 

the right. If indexes 10 through 13 are expressed as 
25 binary numbers, the operation can be implemented 

simply by checking the code of the least significant bit 
of the index on the left. 

In this embodiment, as shown in Equations (5) above. 
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the indexes of channels 1 through 3 are changed according 
to the index of another channel. For example, index jl 
of channel 1 changes according to index iO of channel 
0, so that icil[jl] = {1, 9, 11, 25,} when iO = 0, and 
5 icil[jl] = {5, 13, 21, 29} when iO = 1. 

FIG. 2 and FIG. 3 are flowcharts showing an example 
of a pulse search algorithm for each channel in a. coding 
method according to this embodiment. 

In FIG. 2 and FIG. 3, loop 0 is a loop in which iO 
10 is changed from 0 through 1, loop 1 is a loop in which 

11 is changed from 0 through 3, loop 2 is a loop in which 

12 is changed from 0 through 3, and loop 3 is a loop in 
which 13 is changed from 0 through 3. 

In FIG. 2 and FIG. 3, first, iO, il, and 12 are fixed 
15 at 0, and as the first stage, y and H in each 13 are 
calculated in loop 3, and maximum values ymax and Hmax 
thereamong, and 10, il, 12, and 13 at that time are stored 
as iiO, ill, 112, and 113 respectively. In this case, 
the channel pulse positions searched for are ici3[j3] 
20 = {3, 11, 19, 27} . 

Next, as the second stage, 12 is incremented in loop 
2, and the above first-stage computations are performed 
for each 12. When iO = 0, 11 = 0, and 12 = 1, the channel 
3 pulse positions searched for in the first stage are 
25 ici3[j3] == {7, 15, 23, 31}. Thus, the channel 3 pulse 
positions searched for in the first stage change according 
to the values of iO, 11, and 12. 

Then, as the third stage, 11 is incremented in loop 
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1, and the above first-stage and second-stage 
computations are performed for each il. In this case, 
the channel 2 pulse positions searched for in the second 
stage change according to the values of iO and il. 
5 Lastly, as the fourth stage, iO is incremented in 

loop 0, and the above first-stage, second-stage, and 
third-stage computations are performed for each iO. In 
this case, the channel 1 pulse positions searched for 
in the third stage change according to the value of iO. 

10 Thus, in this embodiment , us ing an n- fold loop sear ch 

algorithm (where n is the number of channels), internal 
loop candidate positions are changed according to 
loop-external codes . 

Then iiO, iil, ii2, and ii3 are found for which y 

15 and H are largest at all pulse positions searched for. 

As a result, iiO is 3 bits and iil, ii2, and ii3 
are 2 bits each, so that pulse position coding can be 
performed in 9 bits, and together with the polarity codes 
of each channel (1 bit x 4 channels), coding can be 

20 performed with a 13-bit code. Therefore, compared with 
the conventional method, the number of bits necessary 
for coding can be reduced, and a lower bit rate can be 
achieved . 

Meanwhile, 8 locations are possible respectively 
25 for indexes jl, j2, and j3 of channels 1 through 3, and 
therefore there are no positions where there is no pulse 
at all in a subframe, variations of excitation vector 
waveforms corresponding to codes (position information) 



15 

can be secured, and sound quality degradation can be 
prevented . 

Thus, according to this embodiment , pulse positions 
of a predetermined channel are associated with pulse 
positions of another channel by changing the 
predetermined channel index in accordance with another 
channel. As a result, a stochastic excitation vector can 
be represented by fewer bits than heretofore, and 
variations can be secured so that there are no positions 
where there is no pulse at all. 

(Embodiment 2) 

In Embodiment 2, a case is described in which the 
pulse positions themselves of a predetermined channel 
15 are changed in accordance with another channel. 

In this embodiment, channel 0 pulse positions 
iciO[iO], channel 1 pulse positions icil[il], channel 
2 pulse positions ici2 [i2] , and channel 3 pulse positions 
ici3[i3] are as shown below. 
20 iciO[iO] = {4, 7, 12, 15, 20, 23, 28, 31} 

icil[il] = {0, 8, 16, 24} 
ici2[i2] = {2, 10, 18, 26} 
ici3[i3] = {5, 13, 21, 29} 

Here, iO (0^10^7) is the index of channel 0, il 
25 (0^11^7) is the index of channel 1, 12 (0^i2<3) is the 
index of channel 2, and 13 (0^13^3) is the index of channel 
3. 

For example, the i0=0 pulse position is {4}, the 
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iO = l pulse position is {7}, and so on; and the il = 0 pulse 
position is {0}, the il = l pulse position is {8}, and so 
on . 

Then channel pulse positions iciO[iO], icil[il], 
5 ici2[i2], and ici3[i3] are adjusted to kO^ kl, k2 , and 
k3 with indexes iO, il, 12, and 13 by means of Equations 
(6) below. 
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. . . Equation ( 6 ) 
InEquations (6), the''%'' symbol denotes an operation 
that finds the remainder when the numeric value on the 
15 left of (index) is divided by the numeric value on 

the right . 

In this embodiment, as shown in Equations (6) above, 
the pulse positions themselves of channels 1 through 3 
are changed according to another channel. As a result, 
20 adjusted pulse positions kO, kl, k2, and k3 of channels 
0 through 3 are as shown below. 

kO = {4, 7, 12, 15, 20, 23, 28, 31} 

kl = {0, 1, 8, 9, 16, 17, 24, 25} 

k2 = {2, 3, 10, 11, 18, 19, 26, 27} 
25 k3 = {5, 6, 13, 14, 21, 22, 29, 30} 

FIG. 4 and FIG. 5 are flowcharts showing an example 
of a pulse search algorithm for each channel in a coding 
method according to this embodiment. 
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In FIG. 4 and FIG. 5, loop 0 is a loop in which iO 
is changed from 0 through 7, loop 1 is a loop in which 

11 is changed from 0 through 3, loop 2 is a loop in which 

12 is changed from 0 through 3, and loop 3 is a loop in 
5 which i3 is changed from 0 through 3. 

In FIG. 4 and FIG. 5, first, iO, il, and 12 are fixed 
at 0, and as the first stage, y and H in each 13 are 
calculated in loop 3, and maximum values ymax and Hmax 
thereamong, and iO, il, 12, and 13 at that time are stored 
10 as iiO, ill, 112, and 113 respectively. 

Next, as the second stage, 12 is incremented in loop 
2, and the above first-stage computations are performed 
for each 12 . 

Then, as the third stage, il is increased in loop 
15 1, and the above first-stage and second-stage 
computations are performed for each 11. 

Lastly, as the fourth stage, 10 is increased in loop 
0, the above first-stage, second-stage, and third-stage 
computations are performed for each 10, and iiO, ill, 
20 112, and ii3 are found for which y and H are largest at 
all pulse positions searched for. 

As a result, iiO is 3 bits and ill, 112, and 113 
are 2 bits each, so that pulse position coding can be 
performed in 9 bits, and together with the polarity codes 
25 of each channel (1 bit x 4 channels), coding can be 

performed with a 13-bit code. Therefore, compared with 
the conventional method, the number of bits necessary 
for coding can be reduced, and a lower bit rate can be 
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achieved . 

Meanwhile, 8 locations are possible respectively 
for the adjusted pulse positions (kl, k2, and k3) of 
channels 1 through 3, and therefore there are no positions 
5 where there is no pulse at all in a subframe, variations 
of excitation vector waveforms corresponding to codes 
(position information) can be secured, and sound quality 
degradation can be prevented. 

Thus, according to this embodiment, by changing the 

10 pulse positions of a predetermined channel in accordance 
with another channel, a stochastic excitation vector can 
be represented by fewer bits than heretofore, and 
variations can be secured so that there are no positions 
where there is no pulse at all. 

15 In a stochastic codebook provided in a speech 

decoding apparatus, a stochastic excitation vector 
searched for by a speech coding apparatus can be found 
by performing computations by means of an above-described 
search algorithm on codes of each channel coded and 

20 transmitted in an above-described embodiment. 

In the above embodiments, a 2's remainder is found 
as variations are assumed to be 2-fold, but the present 
invention is not limited to this, and is also effective 
in a case where the numeric value for which a remainder 

25 is found is made larger, to 3 or more, in order to achieve 
a still lower bit rate and extended subframe length. 

Also, in the above embodiments, information of a 
plurality of channels is integrated by means of addition. 
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but the present invention is not limited to this, and 
is also effective in a case where a more sophisticated 
function, such as weighted addition (addition with 
multiplication byaconstant) orarandom number generator, 
5 is used. 

Furthermore, in the above embodiments, a value 
reflecting information of another channel is extracted 
by means of multiplication, but the present invention 
is not limited to this, and is also effective in a case 

10 where a more sophisticated function is used, such as when 
a random number generator or conversion table is used. 

Moreover, in the above embodiments, a case has been 
described in which an algebraic codebook is used and an 
impulse position corresponds to a code, but the present 

15 invention is not limited to this, and is also effective 
in a case where a stochastic codebook is composed of sums 
of partial waveforms, and the starting position thereof 
corresponds to a code. 

Also, in the above embodiments, a case has been 

20 described in which an algebraic codebook is used and an 
impulse position corresponds to a code, but the present 
invention is not limited to this, and is also effective 
in a case where a stochastic codebook is composed of a 
multiplicity of fixed waveforms stored in ROM, and an 

25 excitation vector waveform is created by the sum of a 
plurality thereof, and that waveform number corresponds 
to a code. In this case, the present invention can be 
applied easily by replacing ''position" with ''waveform 
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number . " 

As is clear from the above description, according 
to the present invention, by performing coding with a 
pulse position of a predetermined channel associated with 
5 a pulse position of another channel, and taking a code 
combining this and a polarity code as a stochastic codebook 
excitation vector code, it is possible to represent a 
stochastic excitation vector with fewer bits than 
heretofore, and to secure variations so that there are 
10 no positions where there is no pulse at all. 

This application is based on Japanese Patent 
Application No . 2 002-3307 68 filed on November 14, 2002, 
the entire content of which is expressly incorporated 
by reference herein. 
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Industrial Appl icabi 1 i t y 

The present invention is applicable to a CELP speech 
coding apparatus /speech decoding apparatus. 
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